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We briefly review the status of the calculation of next-to-next-to-leading order corrections to large angle Bhabha 
scattering in pure QED. In particular, we focus on the analytic calculation of the two- loop virtual corrections 
involving a heavy-flavor fermion loop, which was recently completed. We conclude by assessing the numerical 
impact of these corrections on the Bhabha scattering cross section at colliders operating at a center of mass energy 
of 1 GeV and at the future ILC. 



1. Introduction 

Bhabha scattering [1] is the scattering pro- 
cess of an electron and a positron going into an 
electron-positron final state. The scattering of 
electrically changed particles is always accompa- 
nied by electromagnetic radiation, so that it be- 
comes necessary to consider Bhabha scattering 
events of the type e + e~ — > e + e~ + rvy. 

Since Bhabha scattering has a large cross sec- 
tion and involves charged leptons in the final 
state, it is a process that can be measured ex- 
perimentally with very high accuracy. For ex- 
ample, in the KLOE/DA<f>NE experiment with a 
center of mass energy of yfs ~ 1 GcV, the cross 
section integrated over the allowed phase space 
region with scattering angle 55° < 9 < 125° is of 
about 400 nb, known with 1 permille accuracy. 

Both at high energy e + — e~ colliders and at 
e + — e~ colliders operating at a intermediate cen- 
ter of mass energies (y/s ~ 1 — 10 GeV), it is 
possible to find kinematic regions in which the 
Bhabha scattering cross section is dominated by 
QED. This is the reason why the Bhabha scat- 
tering cross section can be calculated with high 
precision in perturbation theory, including next- 
to-next-to-leading order (NNLO) corrections. 

Due to the nature of the two properties listed 
above, the Bhabha scattering was chosen for the 
luminosity determination at LEP, at the flavor 



*Work supported by the Swiss National Science 
Foundation (SNF) under contract 200020-117602. 
Zurich preprint number ZU-TH 11/08. 



factories (DA<f>NE, VEPP-2M, etc.), and it will 
be employed for the luminosity measurement at 
the future International Linear Collider (ILC). 
At colliders operating at a center of mass energy 
yfs ~ 1 — 10 GeV, the luminosity is measured 
by observing Bhabha events at large scattering 
angles, 9 ~ 90°; in the past, at colliders operat- 
ing at high center of mass energy such as LEP 
(y/s ~ 90,200 GeV), the luminosity was mea- 
sured by observing Bhabha events at small scat- 
tering angles, 9 ~ 1°. At the ILC (y/s ~ 500 
GcV) , in order to disentangle the luminosity spec- 
trum, it will be necessary to control both regions, 
6~ 1° and<9~ 90°. 

The luminosity £ is defined as the ratio be- 
tween the rate of events of a reference scatter- 
ing process divided by the theoretical cross sec- 
tion of the same process: C = N/a t h- It en- 
ters, as a normalization factor, in the measure- 
ments of all other cross sections. Therefore, it 
has to be known with a very high accuracy. Since 
SL = SN+5ath, and SN is to a large extent under 
control, a crucial point is to reduce the theoreti- 
cal uncertainty 5a th to the greatest extent. This 
involves several steps, that we list schematically 
below. 

i) To begin, one needs a precise knowledge of 
the matrix elements, that is, to control the 
perturbative corrections to the basic pro- 
cess e + er — ► e + e~ + ny. At the level of 
accuracy required by the experimental mea- 
surement, the NNLO corrections have to be 
taken into account. Therefore, from a dia- 
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grammatic point of view, it is necessary to 
evaluate two-loop 2 — > 2 Feynman diagrams 
as well as one-loop 2-^3 and tree-level 
2^4. 

ii) Secondly, the theoretical cross section has 
to be inclusive of the soft-photon radiation 
and hard-collinear photon radiation. 

iii) Finally, it is necessary to consider the spe- 
cific experimental set up. Usually detectors 
do not cover the whole solid angle but just 
a portion of it; when calculating the cross 
section, it is necessary to account for these 
experiment specific geometric constraints. 

An accurate theoretical description of the Bhabha 
scattering must consider all these aspects. For 
this purpose, the best tools are Monte Carlo event 
generators (MC) that merge the fixed-order cal- 
culations with parton showers that take into ac- 
count collinear and soft electromagnetic radia- 
tion. In the past, many MCs for the evaluation 
of the Bhabha scattering cross section were de- 
veloped [2 3 4I5I6I7I8] . For a detailed discussion 
of their properties see [9] . 

In the rest of this proceeding, we will concen- 
trate on the analytic calculation of the perturba- 
tive corrections to the matrix element. 

2. Perturbative Corrections 

Since the theoretical error on the Bhabha scat- 
tering differential cross section directly affects 
the precision of the luminosity determination, in 
recent years a significant effort was devoted to 
the calculation of perturbative corrections to this 
scattering process. The NLO corrections are well 
known in the full Standard Model [TO]. For what 
concerns the electroweak corrections at NNLO, 
only the logarithmic enhanced terms are known 
[llj . In pure QED, the situation is significantly 
different. The first complete diagrammatic calcu- 
lation of the two-loop QED virtual corrections to 
Bhabha scattering can be found in |I2j . However, 
this result was obtained by setting the electron 
mass m e to zero from the start, and by employ- 
ing dimensional regularization (DR) to regulate 
both soft and collinear divergencies. 



Today, the complete set of NNLO corrections to 
Bhabha scattering in pure QED have been eval- 
uated using m e as the collinear regulator, as re- 
quired in order to include these fixed order cal- 
culations in MCs. The two-loop Feynman dia- 
grams involved in the calculation can be divided 
in three gauge independent sets: i) diagrams 
without fermion loops ("photonic"), ii) diagrams 
involving a closed electron loop, and iii) diagrams 
involving a closed loop of a fermion different from 
the electron. 

A large part of the NNLO "photonic" correc- 
tions can be obtained in a closed analytic form, 
including the full dependence on the electron 
mass, using a technique that has by now become 
standard in multi-loop calculations. This tech- 
nique is based upon the Laporta algorithm |I3j for 
the reduction of the Feynman diagrams to Master 
Integrals (Mis), and further based upon the dif- 
ferential equation method [14] for the analytical 
evaluation of the latter. With this technique, it 
is possible to calculate the NNLO corrections to 
the form factors [15] (see also the analogous cases 
for heavy quarks jTB] and Higgs [17]) and to pro- 
vide the photonic corrections to the cross section 
with the exception of the ones originating from 
two- loop boxes [18]. The calculation of the two- 
loop photonic boxes retaining the full dependence 
on m e is beyond the reach of this technique (par- 
tial results can be found in [TO]). However, the 
full dependence on m e is not phcnomcnologically 
needed [18]. Fortunately, the physical problem 
exhibits a well defined mass hierarchy. The mass 
of the electron is always very small compared with 
other kinematic invariants, and it can be safely 
neglected everywhere except when it acts as a 
collinear regulator. The collinear structure of the 
NNLO corrections is the following: 




where L e = In (s/rn^.) and where £W are func- 
tions of 9 and of the mass of the heavy fermions 
involved in the virtual corrections. The approx- 
imation given by Eq. ([1]) is sufficient for a phe- 
nomenological description of the process. 

In the case of photonic corrections, the cocffi- 
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cients of the square and single collinear logarithm 
in Eq. ([!} were obtained in |20I21| . The pre- 
cision required for luminosity measurements at 
e + e~ colliders demanded the calculation of the 
non-logarithmic coefficient, that was obtained in 
22J through the infrared matching to the massless 
approximation. The technique of [22j allowed the 
reconstruction of the photonic differential cross 
section in the s 3> m\ ^ limit from the calcu- 
lation in [12] , where m e was set to zero from the 
start. The method employed in [22] involves a 
change of regularization scheme for the collinear 
divergencies originating from a vanishing electron 
mass. A method based on a similar principle was 
subsequently developed in [23124) : the authors of 
[24] confirmed the result of [22] . 

The corrections involving a closed electron loop 
were calculated diagrammatically with the La- 
porta algorithm and the differential equation 
method in [25) . The analytic result, which re- 
tains the full dependence on m e , was expressed 
in terms of one- and two-dimensional Harmonic 
Polylogarithms (HPLs) of maximum weight three 
[26j . By expanding the cross section in the limit 
s,\t\ ^> TOg, the ratio of the NNLO corrections 
to the Born cross section can be written as in 
Eq. (JTJ. However, the series starts with a cubic 
instead of a square collinear logarithm. As ex- 
pected, the cubic collinear logarithm cancels out 
once the contribution of the soft pair production 
graphs is added to the virtual and soft corrections 
[2 7) . The analytic expressions of the expansion of 
the results of [25] was confirmed in |28|24j . 

3. Heavy Flavor Corrections 

We consider now the corrections originating 
from two-loop Feynman diagrams with a heavy 
fermion loop. Since this set of corrections in- 
volves an additional mass scale with respect to 
the corrections analyzed in the previous section, 
a direct diagrammatic calculation is in principle 
a more challenging task. Recently, in [24] the 
heavy flavor NNLO corrections to the Bhabha 
scattering cross section were obtained in the limit 
s, \t\, \u\ >• m 2 ^> m\, where nij is the mass of 
the heavy fermion in the loop. This result was im- 
mediately confirmed in [55] . However, the results 



obtained in this approximation cannot be applied 
to the case in which the i/s < nrif (such as in the 
case of tau loops at yfs ~ I GeV), and they apply 
only to a relatively narrow angular region perpen- 
dicular to the beam direction when y/s is not very 
much larger than mf (such as in the case of top 
quark loops at 1LC). ft was therefore necessary to 
calculate the heavy flavor corrections to Bhabha 
scattering assuming only that s, mj 3> m 2 e . 

The technical problem is simplified by carefully 
considering the structure of the collinear singu- 
larities of this set of corrections. The ratio of the 
NNLO heavy flavor corrections to the Born cross 
section is given by 




It is possible to prove that, in a physical gauge, all 
the collinear singularities factorize and can be ab- 
sorbed in the external field renormalization [29] . 
This observation has two important consequences 
in the case at hand. The first consequence is 
that box diagrams are free of collinear divergen- 
cies in a physical gauge; since the sum of all boxes 
forms a gauge independent block, it can be con- 
cluded that the sum of all box diagrams is free 
of collinear divergencies in any gauge. The sec- 
ond consequence is that the single collinear log- 
arithm in Eq. ^ arises from vertex corrections 
only. Moreover, if one chooses on-shell UV renor- 
malization conditions, the irreducible two loop 
vertex graphs are free of collinear singularities. 
Therefore, among all the two-loop diagrams con- 
tributing to the NNLO heavy flavor corrections 
to Bhabha scattering, only the reducible vertex 
corrections of Fig. [T] are logarithmically divergent 
in the m e — > limild, and they can be easily 
calculated even though they involve two different 
mass scales. By taking advantage of these two 
facts we were recently able to obtain the NNLO 
heavy flavor corrections to the Bhabha scatter- 
ing differential cross section [30], assuming only 
that s, |t|, |ti|,m? ^> m\. In particular, in obtain- 
ing the analytic expression for the NNLO cross 

2 Additional collinear logarithms arise also from the inter- 
ference of one-loop vertex and one-loop self-energy dia- 
grams. 
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Figure 1. Reducible two-loop vertex diagram con- 
tributing to the heavy flavor corrections. These kind 
of diagrams contain single electron mass logarithms. 





Free of 
collinear 
poles 



Free of 
collinear 
poles 



Figure 2. Cancellation of the collinear poles 
among one- and two-loop box diagrams in the 
Feynman gauge. The diagrams are calculated by 
setting m e — from the start. 



section, we worked in the Feynman gauge, set- 
ting m e = from the start in all the diagrams 
with the exception of the reducible ones of the 
kind shown in Fig. Q] and in the interference of 
one-loop graphs. This procedure allowed us to 
effectively eliminate a mass scale from the two- 
loop boxes. The latter could then be evaluated 
with the techniques already employed in the di- 
agrammatic calculation of the electron loop cor- 
rection^]. In this approach, individual box dia- 
grams are singular in the m e — > limit (since 
they are calculated in the Feynman gauge); their 
collinear singularities appear as additional poles 
in the dimensional regulator e. However, it is easy 
to prove that such divergencies cancel in the sum 
of all the box diagrams, as schematically illus- 
trated at the one- and two-loop level in Fig. [2j By 
expanding the analytic results of [30j . it was pos- 
sible to check the heavy flavor cross section in the 
s, \ t\, \u\ S> ml 3> nig limit, which was previously 
known [24 28]. At intermediate energy colliders 
such as DA$NE, the exact dependence on m/ of 
the results of [30] allows us to account for the 



contribution of muons, taus, b- and c-quark loops 
to the Bhabha scattering cross section. At high- 
energy colliders such as the ILC, the top quark 
contribution can also be exactly evaluated. In the 
case in which the heavy flavor fermion is a quark, 
it was simple to modify the calculation of the two- 
loop self-energy diagrams to obtain the mixed 
QED-QCD corrections [30]. An alternative nu- 
merical approach to the calculation of the heavy 
flavor corrections to Bhabha scattering, based on 
dispersion relations, was pursued in [33] , With 
this method it is also possible to evaluate the con- 
tribution of the light quarks vacuum polarization 
to the Bhabha scattering cross section, that can 
be obtained by convoluting the kernel functions 
with the data concerning the cross section of the 
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hadrons. 



3 The necessary Mis can be found in 130131132] 



4. Conclusions 

The numerical impact of the photonic and elec- 
tron loop QED corrections to the Bhabha scat- 
tering cross section at flavor factories was care- 
fully examined in [8.9J, in the context of the MC 
BABAYAGA. A similar analysis of the heavy 
flavor NNLO corrections is not yet available. 
However, it is possible to evaluate numerically 
the NNLO heavy flavor corrections and compare 
them with the other contributions to the cross 
section. This can provide an estimate of the 
numerical impact of the heavy flavor corrections 
once these are introduced in a MC and correctly 
matched with the parton shower. In Table [T] 
(above) we show the results of such an evaluation 
for = 1 GeV and for 50° < 6 < 130° (see [30] 
for details). The muon loop diagrams dominate 
the heavy flavor corrections; they are an order of 
magnitude larger than the corrections involving 
heavier fermions. They reach 1/2 permille of the 
Born cross section at large 8. In the lower part of 
Table [T] we show the results concerning the case 
in which y/s — 500 GeV. In this case, the muon 
and the tau contributions are of the same order, 
while the top contribution is suppressed by an 
order of magnitude at 90° and by four orders of 
magnitude at 1°. 

In conclusion, the calculation of the two-loop 
corrections to Bhabha scattering in QED is now 



> 
o 



> 

a; 

a 

o 
o 

US 



e 


phot (10~ 4 ) 


e (10- 4 ) 


a (io- 4 ) 


c (10- 4 ) 


T (10- 4 ) 


b (io- 4 ) 


50° 


36.688225 


17.341004 


1.7972877 


0.3605297 


0.0264013 


0.0040026 


70° 


41.240039 


19.438718 


2.6504950 


0.5795114 


0.0465329 


0.0065839 


90° 


45.780639 


21.463240 


3.4581845 


0.6528096 


0.0576348 


0.0077240 


110° 


49.366078 


23.099679 


4.0922189 


0.5082196 


0.0495028 


0.0065277 


130° 


50.349342 


23.847394 


4.4392717 


0.2421310 


0.0273145 


0.0039094 
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phot (10~ 3 ) 


e (10~ 3 ) 


/i (io- 3 ) 


t (IO" 3 ) 


t (io- 3 ) 




1° 


7.0074592 


3.4957072 


0.9690710 


0.1542329 


0.0000575 




50° 


14.819671 


7.5740980 


2.3185800 


1.8411736 


0.1707137 




70° 


15.687591 


8.0081541 


2.3708714 


2.0072240 


0.2998535 




90° 


16.560845 


8.4172449 


2.4207950 


2.1521199 


0.4202418 




110° 


17.270026 


8.7451035 


2.4090920 


2.2456055 


0.4979010 




130° 


17.512918 


8.8954702 


2.2543834 


2.2446158 


0.5287459 



Table 1 

The second-order photonic, electron, muon, c-quark, r-lepton, 6-quark and i-quark contributions to the 
differential cross section of Bhabha scattering at \fs = 1 GeV and y/s = 500 GeV in units of 10~ 4 (above) 
and 10~ 3 (below) of the Born cross section. The c, b and t contributions include the 0(aa s ) term. 



complete; these results remove the last piece of 
pure theoretical uncertainty in the luminosity de- 
termination at low- and high-energy accelerators. 
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